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ABSTRACT 


tA  spectroscopic  investigation  of  a  plasma  shock  tube  may  reveal 
the  number  density  and  temperature  of  the  constituents  of  the  plasma. 

Such  an  investigation  was  conducted  on  a  cylindrical  fused  quartz  shock 
tube  and  a  square,  three  inch  aluminum  shock  tube  a  I  ihi  Phygl-pal 
Laboratory-,-  Aarosppce  Corporation,  .Los  Angeles,  -California  in  June-July, 
1961.  Deuterium  and  hydrogen  plasmas  were  investigated  separately.  Gen¬ 
eral  characteristics  including  ringing  frequency  and  longitudinal  visible 
shock  front  velocity  were  determined.  Spectral,  line  analysis  from  a  1.7 
meter  quartz  prism  spectrograph  permitted  the  determination  of  the  shock 
tube  impurities  and  relative  line  intensities.  Use  of  the  latter  instru¬ 
ment  in  conjunction  with  a  tungsten  ribbon  calibration  lamp  and  a  step 
filter  provided  an  absolute  intensity  standard.  For  the  hydrogen  plasma, 

relative  and  absolute  spectral  line  Intensities  were  obtained  which  per- 
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1.  Introduction 

In  recent  years,  particular  attention  has  been  focused  on  spectro¬ 
scopic  diagnostics  of  plasma  devices.  Early  investigations  and  prelim¬ 
inary  theories  were  based  on  the  results  of  flame,  arc,  and  spark  ex¬ 
periments.  Refinement  to  these  theories  and  formulation  of  new  theories 
have  emerged  from  sophistication  of  the  earlier  plasma  devices  to  include 
among  others  that  of  the  plasma  shock  tube.  Within  such  a  shock  tube, 
the  plasma  is  formed  by  transferring  energy  from  a  shock  generator  to 
the  gas  filling.  The  presence  of  the  shock  wave  and  subsequent  energy 
transfer  is  evident  from  the  appearance  of  a  luminous  shock  wave  within 
the  tube,  the  properties  of  this  wave  being  the  subject  of  the  spectro¬ 
scopic  studies. 

The  two  spectroscopic  instruments  used  in  this  shock  tube  investi¬ 
gation  were  the  streak  camera  and  the  spectrograph.  From  these  instru¬ 
ments,  the  derived  experimental  data  permitted  determination  of  the 
following  plasma  parameters: 

a.  shock  tube  impurities 

b.  luminous  shock  front  velocity 

c.  electron  and  ionic  temperatures 

d.  electron  and  ion  number  densities 

Specific  theories  used  to  calculate  these  parameters  Include  those  of 
relative  and  absolute  line  intensities,  line  broadening,  and  apparent 
series  limit.  Throughout  this  paper,  we  have  briefly  mentioned  each 
theory  prior  to  applying  the  pertinent  results  derived.  A  more  complete 
explanation  of  these  theories  will  be  found  in  the  references  provided. 

The  actual  investigation  was  conducted  at  the  Aerospace  Corporation 
Physical  Research  Laboratory  on  two  separate  shock  tubes,  each  briefly 
described  in  the  next  section.  Calculations  were  initiated  at  Aerospace 
Corporation  and  completed  at  the  U.S.  Naval  Postgraduate  School. 

It  should  be  pointed  out  that  not  all  of  the  previously  mentioned 
plasma  parameters  were  determined  for  each  shock  tube.  Our  study  is 
only  intended  as  a  preliminary  investigation  to  future  additional  and 
more  extensive  studies  to  be  performed  on  both  shock  tubt s. 
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2.  Experimental  Set-up 

Components  comprising  the  experimental  set-up  may  be  generally  div¬ 
ided  into  shock  tubes,  magnetic  piston  driver  and  related  power  components, 
gas  supply  and  vacuum  system,  and  monitor  and  test  equipment.  A  photo¬ 
graph,  Fig.  1,  shows  the  physical  placement  of  the  experimental  set-up 
during  a  spectroscopic  investigation.  A  schematic  of  the  experimental 
setup  is  shown  in  Fig.  2.  In  this  section,  only  a  brief  description  will 
be  given  of  the  components  mentioned  above. 

A.  Shock  tubes 

The  first  shock  tube  placed  under  investigation  was  a  General 
Electric,  type  204,  fused  quartz  tube.  Approximate  dimensions 
were  45  inches  in  length  with  inside  and  outside  diameters  of 
74  and  79  mm,  respectively.  A  chemical  analysis  by  the  manu¬ 
facturer  was  conducted  to  determine  tube  material  impurities 
(see  Table  I)  .  Prior  to  making  a  spectroscopic  investigation, 
the  tube  was  modified  by  the  placement  of  two  replaceable 
one-inch  diameter  flat  optical  quality  sapphire  windows  about 
10  and  20  centimeters  from  the  driver  assembly.  The  latter 
modification  prevented  the  necessity  of  making  spectral  ob¬ 
servations  through  cloudy  tube  inner  surface  impurity  film 
deposits . 

The  second  shock  tube  was  rectangular  in  shape  and  was  fabri¬ 
cated  of  aluminum  by  Aerospace  Corporation.  Approximate  di¬ 
mensions  were  an  overall  length  of  36  inches  and  an  inner 
and  outer  dimension  of  three  and  four  Inches  square  respective¬ 
ly.  Spectroscopic  observations  were  made  through  six  inch 
diameter  optical  quality  fused  quartz  windows  12.6  inches  from 
the  driver  assembly.  The  smaller  windows  as  noted  in  Fig.  3 
closer  to  the  driver  assembly  were  for  later  investigations 
(interferometric)  to  be  conducted  on  the  shock  tube. 

B.  Magnetic  Piston  Driver  and  Related  hower  Components 

The  magnetic  piston  driver  was  a  modified  form  of  the  "conical 
shock  tube"  type^.  Power  for  the  magnetic  field  driver  was 
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derived  from  a  40  KV  potential  between  a  capacitor  bank  of 
eight  General  Electric  (type  14P97)  l.O^uf  can  capacitors 
(25  KV)  and  a  floating  ground  strap  at  a  negative  15  KV 
potential  (see  Fig.  2.). 

C.  Gas  Supply  and  Vacuum  System 

The  vacuum  system  consisted  of  a  fore  pump,  diffusion  pump, 
and  liquid  nitrogen  trap,  an  arrangement  providing  pressures 
as  low  as  10  ^  cm  Hg.  Vacuum  gauges  used  were  a  Veeco,  type 
RG-2A  (pressures  less  than  100>4.Hg)  and  a  Autovac,  type  3294B 
(pressures  above  lOO^Hg).  In  addition,  a  purge  pump  was  in¬ 
corporated  into  the  gas  supply  line.  The  gas  supply  system 
consisted  of  bottled  gas  and  associated  tubing  and  valve  com¬ 
ponents. 

D.  Monitor  and  Test  Equipment 

The  major  components  of  moniter  and  test  equipment  were: 

1)  Gaertner  Scientific  Corporation  Quartz  Spectrograph, 
type  L  254 

2)  Philips  Research  Laboratories  Tungsten  Ribbon  Lamp, 
type  W2  (with  intensity  calibration) 

Additional  equipment  used  in  data  reduction  included: 

1)  Leeds  and  Northrup  Recording  Microphotometer,  type 
6/00-2 

2)  McPherson  Instrument  Measuring  Comparator,  type  100 
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3.  Investigation  Procedures 

A  brief  discussion  is  made  below  outlining  the  preliminary  tests  made 
cn  the  shock  tube  cystem.  Mention  is  then  made  of  the  procedures  followed 
in  the  subsequent  spectroscopic  investigation. 

A.  Preliminary 

The  ringing  frequency  of  the  discharge  circuit  was  determined 
by  placing  a  pickup  coil  about  one-half  inch  above  the  return 
straps  of  the  driver  assembly.  Analysis  of  an  oscilloscope 
Polaroid  picture  of  the  coil  signal  revealed  a  discharge 
circuit  ringing  frequency  of  245  KC. 

The  luminous  shock  wave  produced  in  the  fused  quartz  shock 
tube  was  examined  with  a  rotating  mirror  streak  camera. 

Fig.  4  is  an  enlargement  of  a  typical  luminous  shock  wave 
velocity  smear  picture.  The  appearance  of  weak  shock  waves 
after  the  primary  brijht  shock  may  be  attributed  to  the 
ringing  of  the  discharge  circuit  after  firing.  The  heavy 
dark  lines  in  Fig.  4  are  reference  distance  markers  along 
the  shock  tube  axis.  The  entrance  slit  of  the  streak  camera 
was  aligned  along  the  central  shock  tube  axis  in  order  that 
the  35  mm  film  strip  yielded  a  plot  of  light  intensity  distri¬ 
bution  along  the  tube  axis  versus  time. 

The  velocity  of  the  luminous  shock  front  was  determined  as  a 

function  of  initial  Hydrogen  and  Deuterium  Inlet  pressure  so 

/' 

that  the  optimum  inlet  pressure  for  maximum  shock  velocity 

* 

could  be  obtained  .  Fig.  5  and  6  are  curves  of  Initial  in¬ 
let  pressure  versus  shock  velocity  for  Hydrogen  and  Deuterium 
gases  respectively.  The  optimum  pressures  for  maximum  shock 
velocities  were  determined  to  be  80/t  of  Hydrogen  and  125^  of 
Deuterium  yielding  maximum  shock  velocities  of  17.22  cm/^sec 
and  11.65  cm/^/^sec  respectively. 

*  The  method  used  to  calculate  the  velocity  of  the  luminous  shock  front 
is  shown  in  Appendix  1. 
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B.  Spectroscopic 

Prior  to  firing  either  shock  tube,  the  system  was  evacuated 
to  about  10  cm  Hg  pressure.  The  gas  supply  line  was  then 
purged,  after  which  the  gas  to  be  ionized,  hydrogen  or 
deuterium,  was  flushed  through  the  system  at  approximately 
100  microns  pressure  for  periods  in  excess  of  15  minutes. 

The  gas  pressure  was  then  regulated  to  that  pressure  corres¬ 
ponding  to  maximum  luminous  shock  front  velocity  and  the  tube 
fired.  Between  tube  firings,  flushing  occurred  during  which 
time  the  system  returned  to  the  desired  pressure  without  any 
readjustment  of  the  gas  supply  system. 
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4  Relative  and  Absolute  Line  Intensities 

Both  relative  and  absolute  line  intensities  are  utilized  in  the 

spectroscopic  calculation  of  electron  temperature  and  number  density. 

A  tungsten  ribbon  standard  lamp  exposure  on  the  same  spectroscopic  plate 

with  the  experimental  exposures  provides  the  reference  for  absolute  and 

relative  intensity  calculations.  Densitometer  traces  were  made  of  the 

standard  lamp  and  experimental  exposures.  Accepted  data  reduction  pro- 
2  3  4 

cedures  ’  *  were  then  used  in  the  calculation  of  relative  and  absolute 
line  intensities. 

A.  Relative  Line  Intensity 

In  the  calculation  of  relative  intensities,  the  spectro¬ 
scopic  plate  was  calibrated  using  the  standard  lamp  source 
in  order  to  account  for  the  change  in  emulsion  sensitivity 
with  wavelength  and  the  non-linear  dispersion  of  the  spectro¬ 
graph  used.  A  calibation  curve  was  then  plotted,  Fig.  7, 
which  relates  observed  line  intensities  at  all  wavelengths 
to  a  common  intensity  scale  as  follows: 


X1 

h  I 

observed 

*2 

corrected  B  ^2* 

Xa  1 

(4-1) 


trtiere  B(A.)  is  the  correction  factor  which  depends  on 
emulsion  sensitivity  and  spectrograph  dispersion. 
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The 

corrected  relative 

line  Intensities 

were: 

Element 

Multiplet^’ ^ 

Wavelength 

Relative  Intensity 

normalized  to 

6563 

1.00 

HP 

4861 

0.0282 

% 

4340 

0.00905 

h<5 

4102 

0.00370 

3970 

0.00155 

0  II 

3 

3749 

0.0211 

0  II 

36 

4190 

0.00346 

0  II 

2 

4349 

0.0106 

0  II 

5 

4415 

0.00780 

0  II 

1 

4649 

0.0150 

Ca  II 

4 

3179 

0.0455 

Ca  II 

3 

3737 

0.00308 

Ca  II 

1 

3934 

0.0843 

Ca  II 

1 

3968 

0.0575 

Mg  II 

1 

2796 

1.02 

Mg  II 

1 

2803 

0.906 

Mg  II 

3 

2798 

0.262 

Mg  II 

2 

2936 

0.124 

Ti  II 

2 

3235 

0.105 

Ti  II 

2 

3242 

0.0555 

Tl  II 

1 

3349 

0.118 

Na  I 

6 

5688 

0.00628 

Na  I 

1 

5890 

0.360 

Na  I 

1 

5896 

0.310 

Na  I 

5 

6161 

0.0141 

C  II 

6 

3921 

0.00410 

C  II 

4 

4267 

0.0374 
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B.  Absolute  Line  Intensities 

The  absolute  line  intensities  for  the  Balmer  series  lines 
were  calculated  separately  for  1^  through  H^-  using  cali¬ 
brated  standard  lamp  exposures  for  each  of  these  wavelengths 
in  the  general  method  described  by  Harrison,  Lord,  and 
Loofbourow^. 

Several  assumptions  were  made: 

1)  In  order  to  obtain  a  value  for  the  spectroscopic 

system  transmission  factor,  the  following  representa- 

4 

tive  transmission  values  were  assumed  . 

a.  0.90  for  each  lens 

b.  0.85  for  each  prism 

c.  1.00  for  each  mirror 

In  order  to  minimize  the  error  involved  in  this  assumption, 
the  standard  source  was  used  with  the  same  optical  parameters 
that  were  used  in  the  experiment^ . 

2)  The  path  length  from  the  standard  lamp  to  the  entrance 
slit  was  sufficiently  short  to  avoid  appreciable  hori¬ 
zontal  extinction  of  the  standard  source  radiant  energy^. 

3)  To  obtain  the  radiant  energy,  the  solid  angle  subtended 

4 

by  the  collimator  at  the  slit  was  approximated  by  ; 

W  =  V  (4-2) 

F 

where  W  =>  solid  angle  subtended  by  the  collimator 
at  the  slit 

A  =  area  of  the  collimator 

F  =  collimator  focal  length 

4)  The  vertical  variation  of  the  emulsion  was  negligible. 

5)  The  time  of  radiation  was  approximately  5  x  10  ^  seconds. 
This  time  corresponds  to  the  duration  of  the  bright 
region  in  the  smear  photograph,  Fig.  3. 

g 

6)  The  shock  tube  was  optically  thin  . 
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^  Q* 


The  absolute  energies  and  intensities  calculated  per  shock 
tube  firing  were: 


Line 


H<5 


Energy 

Intensity 

(ergs) 

,  -2  -1  -1. 
(ergs-cn  -sec  -ster  ) 

1.353  x  10"A 

0.6950  x  107 

0.386 

0.1980 

0.133 

0.0683 

0.0288 

0.01478 

0.0220 

0.01130 
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5.  Shock  Tube  Impurities 

The  identification  of  shock  tube  impurities  for  both  shock  tubes, 
although  a  tedious  procedure,  was  of  primary  interest  and  produced  a 
useful  result.  Following  the  extension  of  the  spectrograph  calibration 
curve  to  include  the  2000-7000  £  range,  shock  tube  firings  were  conducted 
under  the  conditions  shown  in  Table  II.  In  conjunction  with  the  above 
mentioned  curve  and  comparator  measurements,  line  wavelengths  were  obtained. 
A  typical  3400  X  center  wavelength  spectrum  is  shown  in  Fig.  8.  Element 
and  multiplet  identification  was  then  made  using  tables  prepared  by 
Charlotte  Moore'*'  .  Line  identification  results  are  shown  in  Table  III 
and  are  valid  for  both  the  fused  quartz  shock  tube  and  fused-quartz-window 
aluminum  shock  tube.  This  impurity  spectra  corresponds  to  the  known  con¬ 
stituents  of  the  shock  tube  electrodes  and  walls  and  verifies  that  ablation 
occurs  during  shock  tube  firing.  This  result  will  prove  useful  in  impurity 
identification  in  similar  shock  tube  systems  employing  fused  quartz  sur¬ 
faces  as  spectral  windows. 
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6.  Electron  Number  Densities 

A  knowledge  of  number  densities  is  necessary  in  several  of  the  methods 
for  determining  electron  temperature.  Two  approaches  to  number  density 
calculation  are  those  of  A.  apparent  series  limit  and  B.  the  Stark  broadening 
profile  comparison.  The  formulation  of  a  hydrogen  plasma  clearly  validates 
the  assumption  of  equating  electron  and  ionic  number  densities. 

A.  Apparent  Series  Limit 

As  the  principal  quantum  number,  m,  of  a  Balmer  line  Increases, 
Stark  broadening  of  that  line  increases.  Therefore,  these 

9 

lines  merge  and  form  an  apparent  contlnum.  Inglis  and  Teller 
derived  a  relationship  between  the  last  discernible  line  before 
the  resulting  apparent  contlnum  and  number  density  for  ap¬ 
plication  to  astrophysical  observations.  A  slight  modification 
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to  their  equation  was  made  by  Lochte-Holtgreven  for  use  with 
radiating  plasma  observations.  This  relationship  is: 

Log  N  ■  23.46  -  7.5  Log  m  (6-1) 

max 


where  N  =  n  =  n.=number  density 
C  X 

principal  quantum  number  of  the  last 
discernible  line  within  the  series. 


m 


max 


Letting  ranax  =  7  for  H^.  (see  Table  ixijresults  in  the  values: 

i  i  m17  "3 

n  =  n,  =  1.3  x  10  cm 
e  i 


B.  Stark  Broadening  Profile 

For  particularly  broad  Balmer  line  profiles,  it  may  be  assumed 
linear  Stark  broadening  is  dominant.  Among  the  initial  treat¬ 
ments  of  this  phenomenon  was  that  of  Holtzmark,  who  considered 

o 

only  the  effects  of  perturbing  ions  on  the  emitting  atoms  . 

As  this  theory  did  not  closely  agree  with  experimental  line 
profile  observations,  further  refinement  of  this  theory  was 
formulated  by  Grlem,  Kolb,  and  Shen  and  included  the  per¬ 
turbing  effects  of  both  electrons  and  ions  on  the  hydrogen 

and  are  shown 


atoms.  The  experimental  line  profiles  of  ^ 
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in  Fig.  9  and  10,  each  a  plot  of  normalized  line  intensity, 

I/I  ,  versus  displacement  from  center  wavelength, A  A* 

°  10 
Theoretical  profiles,  as  determined  by  Griem,  Kolb,  and  Shen 

are  also  shown  superimposed  on  the  experimental  profiles.  In 

plotting  the  theoretical  profiles,  the  useful  relationships 

i  j  10 

developed  were  : 

AA  =  1-25  X  109  N2/3  a  (6-2) 


i/l  =  ■(«««)■  8  (0) 


0.625  x  10’9  N2/3 


(6-3) 


where  N  =  electron  number  density. 

C(  =  ratio  of  the  distance  from  the  unperturbed 

line, A  A,  Co  the  Holtzmark  normal  field 

,2/3 

S(a)  =  line  profile  function  such  that  j  S(cr)  da  =  1. 

-o-o 


strength,  Fq  =  2.61  eN 


-oo 


(area)  =  area  under  the  normalized  experimental  line 
profile. 


Reference  to  Fig.  9  and  10  indicates  that  the  effect  of  tem-  ' 
perature  is  small  and  that  of  number  density  is  predominant  in 
determining  the  theoretical  profile.  From  the  curve  fitting 
then,  it  can  be  seen  that 

1016  <  N  <  1017 

where  N  =  n  =  n,  in  centimeters 
e  i 

In  the  calculations  of  electron  temperature  to  follow,  it  will 

be  assumed  that  electron  and  ionic  number  densities  equal  10*7 

-3 

centimeters  .  This  approximation  is  borne  out  by  the  above 
calculations  and  is  considered  accurate  within  an  order  of 
magnitude . 
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7.  Electron  Temperature  Determination 

The  temperature  of  a  pi.isma  may  be  defined  only  if  there  is  statistical 
11  J2 

or  thermal  equilibrium  *  .  Assuming  thermal  equilibrium,  the  population 

distribution  over  the  excited  energy  levels  varies  according  to  Boltzmann's 
law.  The  spectral  intensity  may  therefore  be  related  to  the  temperature. 

The  intensity  of  an  emission  line  corresponding  to  a  transition  from 
r.  to  m  is  given  by: 


I 

run 


b  h™  hi; 


nm 


BA1™  N  hV 

n  nm 


(7-1) 


where  B  =  instrumental  constant 

W™"  =  N  A™  =  number  of  transitions  taking  place  in  light  source 

second 


N 


.nm 


hi/ 


nm 


=  number  of  atoms  in  initial  state  n 


=  Einstein  transition  probability  = 
=  radiated  photon  energy 


number  of  transitions  of  atom 
second 


For  the  case  of  thermal  or  statistical  equilibrium,  the  population 
distribution  may  be  defined  as; 

Ngn 

Nn  "  u  CXP  ("En/k  Te>  (7’2> 


where  N  =  total  density  of  emitting  species  in  volume  for  which  the 
intensity  is  being  measured. 
gn  =  degeneracy  of  initial  state  n  (statistical  weight) 
u  =  partition  function  of  all  statistical  weights 
En  =  energy  of  the  initial  state  n  above  ground  state 
Tg  =  election  temperature 

Also13,  gn  ArUn=  2/3  x  1016  gjfj/A2  (7-3) 

where  J  =  1  or  2 

X=  wavelength  of  the  transition 

fjs  oscillator  strength  determined  for  absorption 

f^=  oscillator  strength  determined  for  emission 
gn=g2=8tatl8tlcal  weight  of  the  initial  state  n 
gj=  statistical  weight  of  the  final  state  m 
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(7-4) 


Therefore, 

r™  =  OlL  <gn£n)  exp  <-En/k  V 
where  D  =  (SS&)  <2/3  x  1016>  (—  ) 

c 

Equation  7  *4  may  be  used  with  the  absolute  intensity  of  a  single  line 
in  order  to  determine  the  electron  temperature,  but  this  method  is  diffi¬ 
cult  due  to  the  uncertainty  in  the  total  density,  N,  of  the  plasma.  In 

calculating  N,  it  is  necessary  to  determine  the  density  of  each  constituent 

q 

present,  requiring  additional  assumptions  Lo  those  already  madt  .  Also, 
there  is  uncertainty  of  reproducing  the  same  constituent  densities  with 
each  tube  firing  due  to  the  vari'ahle  effects  of  ablation  in  the  shock  tube 


Several  methods  were  attempted  in  the  determination  of  electron 
temperature . 

tiethod  I:  Relative  Intensity 

Equation  7-4  may  be  rearranged  as  follows: 

(7-5) 


Log 


’tE'o'n) 


E  log  e 
n 

kT 


Since  D  is  independent  of  the  energy  E  ,  it  may  take  any  value  and 

n  13 

the  left  hand  side  of  7-5  may  be  expressed  in  arbitrary  units  .  in  the 

calculation  of  temperature  using  the  lines  from  an  element,  the  relative 

13  14 

intensities  and  relative  values  of  (g  f  )  may  be  used  '  .  A  plot  of 


log 


versus  E  for  a  number  of  lines  from  an  element  should  therefore  result 

n  1 
in  a  straight  line  with  slope  =  -  Log  e. 
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Method  II:  Relative  Intensity  Ratios 

Equation  7“4  may  also  be  formed  into  a  ratio  o£  two  lines  from  the 
same  element,  giving  the  resultant  equation: 


(7-6) 


Since  ratios  are  being  used  any  arbitrary  consistent  units  may  be 
used  for  I,y  and  (gf )  .  A  plot  of 


log 


11  Vl  *2f2 


*2^1  8lfl 


versus  (E.  -  E„)  for  a  number  of  pairs  of  lines  from  an  element  should 

1  ^  1 
therefore  result  in  a  straight  line  with  slope  =  -  r-~  log  e  . 

kTe 


In  Methods  I  and  II,  the  resultant  graphs,  a  sample  shown  in  Fig.  11, 
not  only  yield  a  value  for  the  electron  temperature,  but  also  acr  as  a 
quasi-check  for  the  assumptions  of  thermal  equilibrium.  If  the  resultant 
points  do  not  fall  on  a  straight  line  and  the  data  has  been  checked  for 
systematic  errors,  then  any  curvature  indicates  a  lack  of  thermal  equili¬ 
brium  in  the  population  of  the  excited  states  of  the  constituents  and  no 

11  13 

temperature  can  be  defined  (this  curving  was  not  observed  in  our  results).  * 


Method  III  Absolute  Intensity 

The  absolute  intensity  of  the  Balmer  lines  in  emission  may  be  used  in 
the  electron  temperature  determination  as  follows:^ 


abs 


?iie2hf»2S  / 


exp 


-E  +E 

n  o 

kT 


0-7) 
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where  e  =  electron  charge 
h  =  Flunk's  constant 
£  =  oscillator  strength 

n  =  principle  quantum  number  of  the  upper  state 

=  neutral  atom  number  density 

1  3  radiating  layer  thickness 

m  =  electron  mass 
e 

^  3  emission  line  wave  length 
3r  s  energy  of  the  state  n 
Eq  3  energy  of  the  ground  state  3  -  13.595  ev 


The  depression  of  ionization  potential  within  the  plasma  is  a  result 
of  the  presence  of  ions  which  are  a  finite  distance  apart.  The  effect 
of  the  ions  is  to  alter  the  Coulombic  field  in  which  the  electrons  move 
and  may  be  expressed  as:^ 


AE  =  7.0  x  10"7  Ne1/3  ev  (7-8) 

where  Ng  =  electron  number  density. 

17  -3 

This  corresponds  to  a  depression  of  0.32  ev  with  =  10  cm  . 

Assuming  that  the  thickness  of  the  radiating  layer  is  approximately 
given  by  the  shock  tube  diameter***  and  substituting  the  electron  number 
density  determined  earlier,  equations  7-7  and  7-8  reduce  to: 

-1.16  x  104(13.28  +  En)  9) 

Ie_  36.77  +in/l.btA3  \ 

1  »2  *  »J 

where  T  is  in  degrees  Kelvin,  N  in  Centimeters)  3,  Ain  centimeters.  I  . 

-2  -1  -1  0  abs 
in  ergs-cm  -sec  -ster  .,  Er  in  electron-volts,  and  f  and  n  are  dimen¬ 
sionless. 

Equation  7-9  may  be  solved  for  any  of  the  Balmer  lir-a  in  emission  if 
the  neutral  atom  number  density  can  be  obtained.  This  latter  quantity  was 
assumed  equal  to  the  electron  number  density,  that  is  50  percent  ionization. 
The  effect  of  this  assumption  is  discussed  in  the  next  section. 
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Results  of  the  electron  temperature  calculations  are: 


Method 


Element 


Temperature 
degree  Kelvin 


8.  Discussion 

The  calculation  of  accurate  temperatures,  power  levels,  and  densities 
in  a  plasma  is  very  important.  Aside  from  the  ability  to  describe  accurately 
the  physical  parameters  of  the  plasma  system,  these  values  enable  the  further 
derivation  of  experimental  transition  probabilities  and  the  calculation  of 
collision  cross-sections  which  to  date  are  not  sufficiently  known. 

In  association  with  Doctor  Richard  M.  Head,  considerable  time  was 
devoted  to  identification  of  the  various  impurities  present  in  the  shock 
tubes^.  Since  the  observed  impurity  spectrum  was  identical  for  the  quartz 
shock  tube  and  for  the  aluminum  shock  tube  with  quartz  windows,  it  is  felt 
that  the  list  of  spectral  lines  shown  in  Table  III  will  be  of  assistance  in 
'  .purity  identification  for  most  quartz  systems  operating  under  similar  con¬ 
ditions.  It  should  be  pointed  out  that  only  those  elements  that  could  con¬ 
tribute  significantly  to  the  intensity  of  a  particular  line  are  included  in 
Table  III.  There  were,  in  addition  to  the  listed  wavelengths,  some  100  lines 
that  were  of  very  low  intensity  which  are  omitted  due  to  uncertainty  in  their 
origin.  Also,  most  of  the  lines  below  2800  £  are  omitted  because  of  limited 
Identification  data  in  this  region. 

Using  the  spectral  lines  of  the  various  impurities  in  the  calculation 
of  temperature,  a  check  may  be  made  on  whether  a  unique  temperature  can  be 
defined  for  the  plasma 

The  accuracy  of  the  calculated  temperature  and  number  densities  in  the 
shock  tube  is  highly  dependent  on  the  relative  and  absolute  line  intensities 
obtained.  In  order  to  determine  correctly  the  relative  or  absolute  Intensity 
of  a  line,  a  number  of  factors  must  be  carefully  considered.  These  include: 

1)  Spectroscopic  techniques  used 

2)  background  continum  correction 

3)  plate-calibration  corrections  for 

a)  dispersion 

b)  evulsion  sensitivity 

c)  standard  source 

4)  self-absorption  considerations 

5)  computational  errors 
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The  first  three  correction  factors  were  determined  using  standard 
spectroscopic  techniques  described  in  Harrison,  Lord,  and  Loofbourow  , 
and  with  the  special  techniques  mentioned  in  appendix  II. 

The  problem  of  self-absorption,  if  present,  can  seriously  affect  the 
measured  intensities  of  the  lines,  and  therefore  lead  to  meaningless  re¬ 
sults.  Severe  self -absorption  is  characterized  by  flat- topped  lines  or 

separation  of  one  line  into  two  lines  by  strong  absorption  in  the  middle 

13 

of  the  spectral  line  .  For  plasma  systems  operating  at  temperatures 

g 

greater  than  approximately  one  ev,  the  problem  of  self-absorption  is  lessened: 

Such  absorption  is  related  to  the  total  intensity  of  emission  from  an  element, 

which  is  temperature  (energy)  dependent,  and  the  total  number  of  potential 

absorbers.  At  temperatures  (energies)  near  the  maximum  intensity  of  emission, 

the  amount  of  absorption  becomes  negligible  since  the  number  of  absorbers 

present  is  much  less  than  the  number  required  to  absorb  all  of  the  emitted 

energy.  For  hydrogen,  this  temperature  corresponds  to  approximately  one  ev. 

As  severe  self-absorption  was  not  observed,  the  system  was  therefore  assumed 

to  be  optically  thin  for  the  Balmer  series.  Future  experimental  work  with 

this  system  should  include  a  check  on  the  validity  of  the  optically  thin 

plasma  assumption  by  placing  a  mirror  on  the  optical  axis  behind  the  source, 

13 

as  recommended  by  Pearce 

Since  neither  radial  distribution  not  time  resolution  studies  were  con¬ 
ducted  in  this  experiment,  the  intensities  determined  permit  only  the  cal- 

14 

culation  of  a  mean  integrated  temperature  for  the  plasma  . 

Taking  into  account  the  above  considerations,  it  is  believed  that  the 
intensities  calculated  are  representative  of  the  populations  of  the  various 
energy  states  and  consequently  may  be  used  in  temperature  and  number  density 
calculations  using  the  Boltzmann  distribution. 

The  two  methods,  apparent  series  limit  and  Stark  Broadening  profile 
comparison,  employed  to  determine  number  densities  are  based  on  the  exist¬ 
ence  of  line  broadening  phenomena.  They  are,  however,  unrelated  as  to  the 
parameters  being  considered. 

In  the  method  of  the  apparent  series  limit,  there  is  no  known  limita¬ 
tion  placed  on  equation  6-1  other  than  a  decrease  in  the  calculated  value 
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by  a  factor  of  two  in  number  density  if  the  plasma  temperature  is  less 

5  8 

than  10  /  ra  (where  m  =  7  in  our  study)  .  This  seems  unlikely  from 
max  max 

the  temperature  calculations  in  this  research. 

In  the  method  of  Stark  broadening  profile  comparison,  the  range  of 

14  -3  19 

applicability  for  experimental  profiles  provided  is  10  £.  n  cm  ^.10 

—  e 

and  10,000<T  °K  ^40,000.  As  previously  mentioned,  the  effect  of  number 

density  variation  is  dominant  over  that  of  temperature  variation.  The 

accuracy  of  the  theoretical  profiles  is  approximately  ten  percent  f  ancj 

the  accuracy  of  experimental  profiles  used  are  estimated  to  be  of  this 

same  order  of  magnitude.  The  experimental  accuracy  is  not  quite  as  high 

for  succeeding  lines  in  the  Balmer  series  because  of  masking  of  the  wings 

17  -3 

by  impurities.  The  choice  of  10  cm  must  be  within  an  order  of  magnitude 

for  electron  number  density  since  it  appears  consistent  with  the  results  of 
17  "3  16  “3  17 

1.3  x  10  cm  and  10  <nfi  cm  <10  by  the  series  limit  and  Stark  broaden¬ 

ing  methods  respectively. 

Prior  to  electron  temperature  calculation  using  absolute  intensities, 

it  was  necessary  to  estimate  the  neutral  atom  number  density.  The  value 
17  -3 

10  cm  was  decided  upon  with  the  stipulation  of  studying  its  variation 
as  a  parameter.  Variation  of  neutral  atom  number  density  by  two  orders 
of  magnitude  affected  the  resulting  electron  temperature  by  less  than  a 
factor  of  two.  In  view  of  this  result,  it  is  not  belived  that  the  elec¬ 
tron  temperature  results  are  materially  affected  by  the  necessity  of  this 
procedure . 

In  the  determination  of  electron  temperature,  the  calculations  were 

limited  by  the  lack  of  transition  probabilities  for  many  of  the  observed 
18 

lines  .  In  methods  I  and  II  the  accuracy  of  the  results  depends  strongly 
on  the  properties  of  the  lines  used  in  the  calculations.  It  is  very  de¬ 
sirable  to  choose  a  number  of  lines  originating  from  a  wide  range  of  upper 

14 

energy  levels  and  occurring  over  a  small  wavelength  interval  .  This  pro¬ 
cedure  was  somewhat  compromised  due  to  the  lack  of  transition  probabilities. 
The  range  of  validity  of  methods  I  and  II  is  up  to  approximately  5  x  10  °K 

g 

with  the  proper  selection  of  the  elements  .  For  each  wavelength  line  there 
is  a  relatively  flat  maximum  in  the  temperature-intensity  profile  which  ex¬ 
tends  over  a  few  thousand  degrees.  Therefore,  when  measuring  a  temperature, 


20 


it  is  more  accurate  to  select  a  line  where  the  slope  of  the  temperature- 
intensity  profile  is  rapidly  changing.  From  data  presented  by  Lochte- 

g 

Holtgreven  ,  the  temperatures  determined  using  Oxygen  II  and  the  Balmer 
lines  were  on  the  rapidly  rising  slopes  of  the  temperature-intensity  pro¬ 
files  and  thus  in  the  more  accurate  range.  Similar  data  were  not  available 
for  the  other  elements. 

17  -3 

Since  the  electron  density  is  relatively  high  (10  cm  ),  there  is  a 

rapid  approach  to  local  thermodynamic  equilibrium.  The  formula  for  re- 

19 

laxation  time  is  given  by  Spitzer  : 


11. A  ,t3/2 

t  <  (T  )  sec. 

N  Z  In  A 


v 


where  In  A  —  20 

A  =  1  for  a  hydrogen  plasma 

Z  a  1  for  a  hydrogen  plasma  . 

„  in17  -3 
N  «  10  cm 

T  =  6270°  K  (0.54  ev) 

Therefore  t  — 10  ^  seconds  or  less,  and  the  concept  of  electron  temperature 

has  some  meaning  in  the  usual  statistical  sense.  This  is  further  verified 

by  the  results  of  calculations  used  in  methods  I  and  II.  In  these  methods 

a  marked  departure  from  equilibrium  is  characterized  by  a  curving,  or  trend, 

13 

in  the  graphic  data  .  This  curving,  and  hence  departure  from  equilibrium, 
was  not  observed  in  our  result  (see  Fig.  11).  Several  other  advantages  were 
gained  by  using  methods  1  and  II.  A  greater  number  of  lines  could  be  used 
with  these  methods  than  in  the  absolute  intensity  method;  and  therefore,  any 
anomally  caused  by  an  Improper  f-value,  self-absorption,  or  intensity  error 
(which  would  result  in  a  wide  variation  from  the  other  valuec)  could  be  easily 
distinguished.  In  addition  the  more  easily  measured  relative  intensity  may 
be  used  instead  of  absolute  intensity,  and  it  is  not  necessary  to  know  ac¬ 
curately  the  composition  of  the  plasma. 
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The  temperatures  calculated  by  methods  I  through  III  varied  from 
0.54  to  4.2  ev.  This  result  is  to  be  expected,  since  each  plasma  con¬ 
stituent  is  likely  to  exhibit  a  different  temperature  during  the  time 
integrated  observation.  Future  investigations  should  include  a  tempera¬ 
ture-time  variation  study  within  the  shock  tubes. 

The  results  of  method  III  using  the  absolute  intensity  agree  closely 
with  the  electron  temperature  calculated  using  the  relative  intensities  of 
the  Balmer  lines.  The  difference  between  the  methods,  which  is  less  than 
a  factor  of  25  percent,  shows  that  the  assumptions  made  in  method  III  and 
those  used  in  the  absolute  intensity  calculations  were  reasonable. 

The  ion  temperature  may  be  estimated  using  the  results  of  the  absolute 

intensity  measurements.  This  problem  has  been  investigated  theoretically 

20 

for  a  pure  hydrogen  plasma  considering  the  following  factors  : 

1)  Energy  transfer  from  atoms  to  electrons  caused  by  elastic 
collisions  between  atoms  and  electrons  and  ions  and  electrons. 

2)  The  energy  dissipation  by  the  electrons  that  is  determined  by: 

a)  inelastic  collisions  of  electrons  with  atoms  and  electrons. 

b)  the  excess  of  emission  over  absorption  in  free-free  transi¬ 
tions  . 

c)  the  excess  of  emission  over  absorption  in  free-bound  and 
bound-free  transitions. 

The  difference  between  electron  and  kinetic  temperatures  under  these  con¬ 
ditions  would  be: 

-3 

T.-T  4  10  T 
i  e  e 

10  -3 

when  N  >  10  cm 
e 

and  5000°K<T  <  105°K 
e 

From  this  relation  the  ion  temperature  is  approximately  equal  to  the 

electron  temperature  for  the  conditions  existing  in  the  shock  tube  (i.e., 

N  =>  10*7  gj,,  ^  an{j  j  >  0.54  ev.). 
e  e  - 
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Conclusions 

A.  Shock  Tube  Characteristics 

Spectroscopic  diagnostic  techniques  reveal  the  following  plasma 
parameters  existing  within  the  two  shock  tubes  with  hydrogen 
Introduced; 

17  -3 

1)  An  electron  and  ionic  number  density  of  approximately  10  cm  , 

2)  A  peak  temperature  of  approximately  4  ev. 

3)  A  condition  of  local  thermodynamic  equilibrium. 

B.  General 

1)  Shock  tube  systems  employing  quartz  windows  and  operating  under 
similar  conditions  are  characterized  by  nearly  identical  im¬ 
purity  spectra. 

2)  Spectral  line  absolute  intensity  measurements  provide  a  useful 
mechanism  for  determining  plasma  temperature. 

3)  Further  investigation  of  the  shock  tubes  under  consideration 
should  Include  time  resolved  and  radial  temperature  and  number 
density  studies. 


Fig.  1.  Experimental  Set-up  During  Spectroscopic  Investigation 


Fig.  2.  Schematic  of  the  Experimental  Set-up 
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Fig.  4.  Shock  Wave  Smear  Photograph 


Wave  Velocity,  cm^Asec 


Fig.  5.  Hydrogen  Pressure  Variation  Effect  on  Shock  Wave  Velocity 
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F  g.  6.  Deuterium  Pressure  Variation  Effect  on  Shock  Wave  Velocity 
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Fig.  7.  Spectroscopic  Plate  Calibration 


Note:  a-plate  1  of  7/18,  b-plate  1  of  7/10;  see  Table  II  for  firing  conditions;  exposures  from  left 
to  right  -  single  shot,  two  shots  superimposed,  -four  shots  superimposed. 


Experimental  and  Theoretical  Stark  Broadening  Profile# 
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Fig.  10.  Experimental  and  Theoretical  1L.  Stark  Broadening 
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Fig.  11.  Electron  Temperature  by  Relative  Intensities  of  Balmer  Series 
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Impurity 

Part 8  per  Million 
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/ 

B 

0.94 

Table  I  -  Fused  Quartz  Tube  Impurities  as 
Provided  by  Manufacturer 
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Table  Ill  Spectral  Line  Identification 


Wavelength 

(X) 

Element 

Multiplet^*^ 

Number 

Excitation  Potential 
Low  High 

Relative 

Intensity 

(visurl) 

2367.064 

A1 

I 

4 

0.01 

5.21 

10 

2395.06 

B 

II 

4 

12.64 

17.79 

6 

2435.160 

Si 

I 

45 

0.78 

5.85 

15 

2478.556 

C 

I 

61 

2.67 

7.65 

10 

2496.773 

B 

I 

1 

0.00 

4.94 

18 

2497.724 

B 

I 

1 

0.00 

4.94 

20 

25l6.896 

Si 

I 

1 

0.01 

4.93 

18 

2514.315 

Si 

I 

1 

0.00 

4.91 

16 

2516.109 

Si 

I 

1 

0.03 

4.93 

20 

2519.203 

Si 

I 

1 

0.01 

4.91 

17 

2524.108 

Si 

I 

1 

0.01 

4.90 

18 

2528.510 

Si 

I 

1 

0.03 

4.91 

18 

2532.38 

Si 

I 

86 

1.90 

6.77 

8 

2541.83 

Si 

III 

8 

10.23 

15.09 

18 

2631.28 

Si 

I 

83 

1.90 

6.59 

12 

2790.768 

Mg 

II 

3 

4.40 

6.83 

4 

2795.523 

Mg 

II 

1 

0.00 

4.41 

22 

2797.989 

Mg 

II 

3 

4.41 

8.83 

9 

2802.698 

Mg 

II 

1 

0.00 

4.40 

21 

2816.189 

A1 

II 

7 

7.39 

11.77 

10 

2852.120 

Mg 

I 

1 

0.00 

4.33 

2 

2881.595 

Si 

I 

43 

0.78 

5.06 

24 

2904.29 

Si 

II 

17 

9.79 

14.04 

6 

2905.70 

Si 

II 

17 

9.80 

14.04 

10 

2936.96 

Mg 

II 

2 

4.41 

8.62 

1 

2987.65 

Si 

I 

1 

0.78 

4.91 

16 

3072.971 

Ti 

II 

5 

0.00 

4.02 

0 

3075.225 

Ti 

11 

5 

0.01 

4.02 

0 

3078.645 

Ti 

II 

5 

0.03 

4.04 

1 

3082.159 

Al 

I 

3 

0.00 

4.00 

12 

3086.225 

Si 

III 

1 

17.63 

21.63 

9 

3088.027 

Ti 

II 

5 

0.05 

4.05 

3 
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Wavelength 

Element 

Multiplet^*^ 

Excitat 

ion  Potential 

Relative 

(X) 

Number 

Low 

High 

Intensity 

(visual) 

3092.716 

A1 

I 

3 

0.01 

4.00 

17 

3092.843 

A1 

I 

3 

0.01 

4.00 

3096.786 

SI 

III 

1 

17.64 

21.62 

1 

3129.76 

Zr 

II 

5 

0.04 

3.98 

O 

3129.16 

Zr 

II 

23 

0.52 

4.47 

4 

3129.44 

0 

II 

14 

25.53 

29.47 

3149.56 

SI 

IV 

2 

26.95 

30.86 

6 

3165.72 

SI 

IV 

2 

26.97 

30.86 

3165.98 

Zr 

II 

5 

0.16 

4.06 

12 

3166.29 

Zr 

II 

48 

0.80 

4.70 

3179.34 

B 

II 

Unci. 

1 

3179.332 

Ca 

II 

4 

3.14 

7.02 

A 

3203.89 

SI 

II 

7 

10.02 

13.87 

6 

3210.04 

SI 

II 

7 

10.03 

13.87 

12 

3210.52 

SI 

III 

Unci. 

3234.517 

Ti 

II 

2 

0.05 

3.86 

A 

3234.00 

Si 

III 

6 

21.62 

25.44 

0 

3236.573 

Tl 

II 

2 

0.03 

3.84 

4 

3239.037 

Ti 

II 

2 

0.01 

3.82 

2 

3241.984 

Ti 

II 

2 

0.00 

3.81 

O 

3241.67 

Si 

III 

6 

21.63 

25.44 

J 

3247.540 

Cu 

I 

1 

0.00 

3.80 

20 

3261.596 

Ti 

II 

66 

1.23 

5.01 

o 

3261.596 

Ti 

II 

89 

1.88 

5.67 

4 

3265.46 

0 

III 

8 

36.32 

40.10 

l 

3273.957 

Cu 

I 

1 

0.00 

3.77 

19 

3279.26 

Zr 

II 

3 

0.09 

3.86 

4 

3282.329 

Ti 

II 

66 

1.22 

4.98 

1 

3284.72 

Zr 

II 

4 

0.00 

3.76 

1 

3302.34 

Na 

I 

2 

0.00 

3.74 

£ 

3302.94 

Na 

I 

2 

0.00 

3.74 

0 

3322.936 

Ti 

II 

7 

0.15 

3.86 

A 

3323.61 

B 

II 

Unci. 

4 

3341.875 

Ti 

11 

16 

0.57 

4.26 

6 

3343.770 

Ti 

II 

7 

0.15 

3.84 

2 

3349.399 

Ti 

II 

1 

0.05 

3.73 

10 

3361.213 

Ti 

II 

1 

0.03 

3.70 

7 

38 


Wavelength 

Element 

Multlplet5’6 

Excitation  Potent  '  )1 

Relative 

<*) 

Number 

Low 

High 

Intensity 

(visual) 

3372.800 

Ti 

II 

1 

0.01 

3.67 

5 

3383.761 

Ti 

II 

1 

0.00 

3.65 

4 

3388.29 

Zr 

II 

2 

0.00 

3.64 

3388.87 

Zr 

II 

74 

0.97 

4.61 

J 

3391.96 

Zr 

II 

1 

0.16 

3.80 

9 

3410.26 

Zr 

II 

11 

0.41 

4.03 

1 

3430.53 

Zr 

II 

11 

0.46 

4.06 

1 

3438.23 

Zr 

II 

1 

0.09 

3.68 

8 

3451.41 

B 

II 

1 

9.06 

12.64 

24 

3452.28 

B 

II 

Unci. 

3479.39 

Zr 

II 

46 

0.71 

4.26 

2 

3481.14 

Zr 

II 

46 

o 

CO 

• 

-> 

4.34 

3 

3483.54 

Zi 

II 

33 

0.75 

4.30 

1 

3496.18 

Zr 

II 

1 

0.04 

3.57 

7 

3505.67 

Zr 

II 

1 

0.16 

3.68 

3505.47 

Zr 

II 

90 

1.53 

5.05 

3 

3504.890 

TI 

II 

88 

1.88 

5.40 

3510.840 

TI 

II 

88 

1.88 

5.40 

1 

3535. 40C 

TI 

II 

98 

2.05 

5.54 

1 

3542.65 

Zr 

II 

113 

1.75 

5.24 

1 

3551.94 

Zr 

II 

1 

0.09 

3.57 

2 

3556.61 

Zr 

II 

9 

0.46 

3.93 

3 

3572.47 

Zr 

II 

1 

0.00 

3.45 

5 

3576.88 

Zr 

II 

9 

0.41 

3.86 

2 

3586.9 

A1 

II 

7  (9)* 

11.80 

15.24 

15 

3590.46 

SI 

III 

7 

21.79 

25.22 

4 

3596.048 

Ti 

II 

15 

0.60 

4.04 

1 

3601.623 

A1 

III 

1 

14.31 

17.74 

5 

3612.352 

A1 

III 

1 

14.31 

17.73 

3 

3614.79 

Zr 

II 

9 

0.36 

3.77 

1 

3685.192 

Ti 

II 

14 

0.60 

3.95 

A 

3685.192 

TI 

II 

14 

0.57 

3.92 

H 

3698.17 

Zr 

II 

71 

1.01 

4.34 

2 

3713.103 

A1 

III 

4 

17.74 

21.07 

1 

*  9  lines  contribute  to  this  multlplet. 
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Wavelength  Element  Multiplet^*^ 
(X)  Number 


3727.33 

0 

II 

3 

3727.72 

Zr 

II 

112 

3736.901 

Ca 

II 

3 

3749.49 

0 

II 

3 

3751.60 

Zr 

II 

71 

3759.291 

Ti 

II 

13 

3759.87 

0 

III 

2 

3761.320 

Ti 

II 

13 

3791.41 

Si 

III 

5 

3791.26 

0 

III 

2 

3796.11 

Si 

III 

5 

3806.56 

Si 

III 

5 

3836.76 

Zr 

II 

16 

3853.657 

Si 

II 

1 

3856.021 

Si 

II 

1 

3862.592 

Si 

II 

1 

3881.97 

Zr 

II 

134 

3882.197 

0 

II 

12 

3905.527 

Si 

I 

3 

3911.960 

0 

II 

17 

3916.997 

C 

II 

4 

3920.677 

C 

II 

4 

3924.44 

Si 

III 

Unci 

3933.664 

Ca 

II 

1 

3944.009 

A1 

I 

1 

3961.523 

A1 

I 

1 

3968.470 

Ca 

II 

1 

3970.074 

\ 

i 

3973.263 

0 

II 

6 

3991.14 

Zr 

II 

30 

3998.98 

Zr 

II 

16 

3998.635 

Ti 

I 

12 

4062.90 

0 

II 

50 

4069.897 

0 

II 

10 

4069.636 

0 

II 

10 

Excitation  Potential 
Low  High 

Relative 

Intensity 

(visual) 

22.88 

1.74 

26.19 

5.05 

1 

3.14 

6.44 

1 

22.90 

26.19 

2 

0.97 

4.26 

1 

0.60 

33.04 

3.89 

36.32 

4 

0.57 

3.85 

3 

21.62 

33.04 

24.88 

36.29 

2 

21.62 

24.88 

4 

21.63 

24.88 

7 

0.56 

3.77 

2 

6.83 

10.03 

9 

6.83 

10.03 

24 

6.83 

10.02 

21 

2.40 

25.55 

5.58 

28.73 

1 

1.90 

5.06 

22 

25.55 

28.71 

1 

16.26 

19.41 

1 

16.26 

19.41 

1 

1 

0.00 

3.14 

18 

0.00 

3.13 

14 

0.01 

3.13 

16 

0.00 

3.11 

14 

10.15 

13.26 

1W 

23.34 

26.45 

1 

0.75 

3.85 

2 

0.56 

0.05 

3.64 

3.13 

1 

28.58 

31.62 

1 

25.53 

25.52 

28.56 

28.56 

3 

40 


Wavelength 

Element 

Multiplet 

<*> 

4072.164 

0 

II 

Number 

10 

4076.00 

C 

II 

36 

4075.868 

0 

II 

10 

4088.863 

SI 

IV 

1 

4097.260 

0 

II 

20,48 

4096.543 

0 

XI 

21 

4101.737 

16 

Si 

1 

4102.926 

I 

2 

4116.104 

Si 

IV 

1 

4121.95 

B 

II 

2 

4128.053 

Si 

II 

3 

4130.884 

Si 

II 

3 

4149.22 

Zr 

II 

41 

4149.897 

A1 

III 

5 

4150.138 

A1 

III 

5 

4189.788 

0 

II 

36 

4267.27 

C 

II 

6 

4267.02 

c 

II 

6 

4275.52 

0 

II 

67 

4340.468 

1 

4349.426 

0 

II 

2 

4347.425 

0 

II 

16 

4366.896 

0 

II 

2 

4414.909 

0 

II 

5 

4416.975 

0 

II 

5 

4479.966 

A1 

III 

8 

4479.891 

A1 

III 

8 

4512.535 

A1 

III 

3 

4529.176 

Al 

III 

3 

4552.654 

Si 

III 

2 

4567.872 

Si 

III 

2 

4641.811 

.  0 

II 

1 

4649 . 139 

0 

II 

1 

4661.635 

0 

II 

1 

4861.332 

1 

Excitation  Potential 
Low  High 

Relative 

Intensity 

(visual) 

25.54 

28.57 

3 

24.27 

25.55 

27.30 

28.58 

6 

23.95 

26.97 

8 

25.72 

25.73 

28.73 

28.74 

2 

10.15 

13.16 

3W 

1.90 

4.91 

2 

23.95 

26.95 

3 

18.60 

21.59 

5 

9.79 

12.78 

23 

9.80 

12.78 

24 

0.80 

3.77 

2 

20.47 

20.47 

23.44 

23. 4^ 

2 

28.24 

31.18 

2 

17.97 

17.97 

20.86 

20.86 

6 

28.73 

31.62 

i 

10.15 

13.00 

10W 

22.90 

25.55 

25.74 

28.39 

2 

22.90 

25.73 

2 

23.34 

23.32 

26.14 

26.11 

3 

20.69 

20.69 

23.45 

23.45 

2 

17.73 

20.47 

1 

17.74 

20.47 

2 

18.92 

21.63 

4 

18.92 

21.62 

•) 

«• 

22.88 

25.54 

1 

22.90 

25.55 

3 

22.88 

25.53 

1 

10.15 

12.69 

15 

41 


Wavelength 

Element 

Multiplet 

(2) 

5041.063 

Si 

11 

Number 

5 

5056.020 

SI 

II 

5 

5105.541 

Cu 

I 

2 

5688.205 

Na 

I 

6 

5682.633 

Na 

I 

6 

5889.953 

Na 

I 

1 

5895.923 

Na 

I 

1 

5957.612 

Si 

II 

4 

5978.970 

Si 

II 

4 

6154.225 

Na 

I 

5 

6160.747 

Na 

I 

5 

6243.36 

Al 

II 

10 

6347.091 

Si 

II 

2 

6371.359 

Si 

II 

2 

6562.817 

Ha 

1 

Excitation  Potential 

Relative 

Low 

High 

Intensity 

(visual) 

10.02 

12.47 

15 

10.03 

12.47 

20 

1.38 

3.80 

1 

2.10 

4.27 

2.09 

4.27 

1 

0.00 

2.10 

18 

0.00 

2.09 

16 

10.02 

12.09 

4 

10.03 

12.09 

10 

2.09 

4.10 

* 

2.10 

4.10 

J 

13.02 

15.00 

1 

8.09 

10.03 

24 

8.09 

10.02 

22 

10.15 

12.04 

25 
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APPENDIX  I 


Streak  Photography 

The  general  arrangement  of  the  streak  camera  Is  shown  In  Fig.  1-1. 


The  object,  0,  to  be  photographed  is 
reflected  off  the  rotating  mirror,  M, 
on  to  the  35  mm  film  strip,  F. 

The  developed  film  strip  pres 
entation  is  shown  in  Fig.  1-2. 

The  velocity  of  the  exposure  on 
the  film  strip  may  be  expressed  as: 
vf  =  2nfR  (1-1) 

where  f  a  frequency  of 

rotating  mirror 

R  -  radius  of  curva¬ 
ture  of  the  streak 
camera  film  holder. 


Fig.  1-1 


Fig.  1-2 


If  the  velocity  of  the  luminous  shock  front,  v,  is  desired  at  point 
A  shown  in  Fig.  1-2,  a  tangent  (i.e.  slope  or  velocity)  at  that  point  may 
be  constructed  as  shown  in  Fig.  1-3.  ^ 

The  right  triangle  shown  has  sides  / 1 

v,  and  mv  where  m  is  the  magnifica-  /  i  / 

f  /  / 
tion  of  the  streak  photograph  and  v  ^  |  S 

is  the  desired  velocity  of  the  lum-  ^  vt 

inous  shock  front  at  point  A.  There-  /  /  I  ®v 


fore: 


-  v,  tan  0 
m  f 


(1-2) 


Combining  1-1  and  1-2,  r  ♦  jJ 

v  -  tan  9  (1-3)  f 

“  Fig.  1-3 

R  was  determined  by  direct  measurement  on  the  streak  camera.  The 
quantities  9  and  m  were  determined  from  the  film  strip,  whereby  f  was  obx* 
tained  through  the  use  of  a  counter  and  timer  connected  to  the  axis  of  the 
rotating  mirror. 
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APPENDIX  II 


Spectroscopic  Techniques 

It  was  necessary  to  make  a  number  of  calibrations  for  the  experimental 
equipment.  Since  the  Gaert.ier  L  254  spectrograph  and  associated  equipment 
will  be  available  at  the  Postgraduate  School  for  further  research  work,  the 
results  of  the  calibrations  made  are  listed  below  for  future  reference. 

A.  Slit  width  calibration 

The  jaws  of  the  slit  were  checked  and  found  to  be  parallel  and  in 
excellent  condlton.  The  various  slit  width  settings  were  measured 
and  plotted.  Fig.  II-l  is  a  curve  of  slit  width  setting  versus 
the  measured  slit  width. 

B.  Timer  Calibration 

The  following  timer  settings  were  determined: 


Dial  Setting 

Measured  Time 

(sec) 

(sec) 

0.5 

0.5 

1.0 

1.0 

1.6 

•  1.5 

2.0 

• 

2.4 

4.0 

4.9 

C.  Step  Filter  Calibration 

The  source  of  the  transmission  values  for  the  step  filter  was 
unknown.  Measurements  made  showed  considerable  variation  with 
the  rated  transmission  values  provided. 


Step  number,  n 

Rated 

VXo 

Measured 

VXo 

1 

1.0 

1.0 

2 

0.627 

0.406 

3 

0.390 

0.244 

4 

0.250 

0.147 

5 

0.156 

0.065 

6 

0.100 

0.030 

7 

0.064 

0.012 
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D.  Other  Calibrations 

The  spectrograph  was  cleaned,  aligned,  and  focused  by  Aerospace 
Corporation  after  receipt  of  the  instrument  from  NOTS,  China  Lake, 
California.  The  Standard  Lamp  and  associated  equipment  was  simi¬ 
larly  calibrated  by  Aerospace  Corporation  prior  to  experimental 
measurements. 

The  measurement  of  radiant  energy  spectroscopically  is  a  very  difficult 
and  sensitive  problem.  There  are  many  factors  involved,  which,  if  not  pro¬ 
perly  compensated  for,  will  lead  to  meaningless  values.  In  general,  the 

2  3  A 

recommended  methods  of  Sawyer  ;  Brode  ;  and  Harrison,  Lord  and  Loofbourow 

were  followed  in  the  intensity  determinations. 

Particular  attention  was  given  to  the  following  items: 

1.  Pre -measurement  checks 

2.  Plate  exposure 

3.  Plate  developing 

A.  Plate  calibration. 

Pre-measurement  checks  included  careful  warming  and  insertion  of  the 
spectroscopic  plate,  check  of  spectrograph  alignment,  focus  and  optical 
set-up,  and  the  complete  darkening  of  the  room  prior  to  measurement. 

Plate  exposures  of  the  shock  tube  firing  and  the  Standard  Lamp  were 
taken  through  identical  optical  set-ups  with  the  exception  of  the  Standard 
Lamp  window  which  necessitated  a  separate  compensation  in  the  data  reduction 
calculations.  The  exposures  were  taken  over  a  short  time  interval  to  reduce 
the  affect  of  ambient  temperature  and  pressure  variations  on  the  spectrograph. 

Temperature  and  time  criteria  were  closely  adhered  to  in  the  developing 
and  processing  of  the  plates.  The  chemicals  used  were  freshly  made,  and 
constant  agitation  was  also  employed. 

The  densitometer  used  in  plate  calibrations  was  carefully  checked  to 
Insure  that  it  was  fully  warmed  and  steadied  out  prior  to  use.  After  each 
densitometer  run  was  made,  a  check  was  immediately  made  to  insure  that  drift 
has  not  occurred  during  the  run. 
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